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Abstract 

The ALICE collaboration reports the measurement of the inclusive }/y/ yield as a function of charged 
particle pseudorapidity density dA^ch/dJ? in pp collisions at y(j = 7 TeV at the LHC. }/y/ particles are 
detected for pt > 0, in the rapidity interval \y\ < 0.9 via decay into e+e^, and in the interval 2.5 < 
y < 4.0 via decay into jJ-^fi ^ pairs. An approximately linear increase of the J/y/ yields normalized to 
their event average (d^Vj/^/dy) / (dA'j/^/dy) with {dNch/drj)/ {dNch/drj) is observed in both rapidity 
ranges, where dNch/dt] is measured within j?]! < 1 and pi > 0. In the highest multiplicity interval 
with (d;V,,h/dTj(bin)) = 24.1, corresponding to four times the minimum bias multiplicity density, 
an enhancement relative to the minimum bias J/\j/ yield by a factor of about 5 at 2.5 < y < 4 (8 at 
\y\ < 0.9) is observed. 
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Understanding the production mechanism of quarkonium states in hadronic collisions is still a challenge 
due to its sensitivity to perturbative and non-perturbative aspects of Quantum Chromodynamics (QCD). 
While the primary production of heavy quark anti-quark (qq) pairs is generally treated as a hard process 
which can be calculated within perturbative QCD, the subsequent formation of a bound colorless qq pair 
is inherently non-perturbative and difficult to treat. The models developed to describe quarkonium pro- 
duction in high energy hadron collisions consequently follow various approaches, mainly differing in the 
relative contribution of the intermediate color singlet and color octet qq states lHHH. Recent theoretical 
work tries to describe consistently E IH [5l the measured production cross section and polarization, in 
particular in light of recent measurements at the LHC ll6ir7l[8ll9l [T0l[m . 

It is also important to consider that a high energy proton-proton collision can have a substantial con- 
tribution from Multi-Parton Interactions (MPI) [12. ,131 . In this case several interactions on the parton 
level can occur in a single pp collision, which can introduce a dependence of particle production on the 
total event multiplicity lfT4l[T5l[T6l . If MPI were mainly affecting processes involving only light quarks 
and gluons, as implemented e.g. in PYTHIA 6.4, processes like J/i// and open heavy flavour production 
should not be influenced and their rates are expected to be independent of the overall event multiplicity. 
However, at the high center-of-mass energies reached at the LHC, there might be a substantial contribu- 
tion of MPI on a harder scale which can also induce a correlation between the yield of quarkonia and the 
total charged particle multiplicity fTT|. An early study that relates open charm production and underlying 
event properties was performed by the NA27 experiment for pp collisions at ^/s = 27 GeV, with the result 
that charged particle multiplicity distributions in events with open charm production have a mean that is 
higher by ~ 20% than the ones without ITSl . 

In |[T9l l20l it has been argued that, due to the spatial distribution of partons in the transverse plane (as 
described in generalized parton distributions), the density of partons in pp collisions will be strongly 
impact parameter dependent. Therefore, the probability for MPI to occur will increase towards smaller 
impact parameters. This effect might be further enhanced by quantum-mechanical fluctuations of the 
small Bj0rken-A; gluon densities. 

The charged particle multiplicities measured in high-multiplicity pp collisions at LHC energies reach 
values that are of the same order as those measured in heavy-ion collisions at lower energies (e.g. they 
are well above the ones observed at RHIC for peripheral Cu-Cu collisions at .^s^ = 200 GeV [21]). 
Therefore, it is a valid question whether pp collisions also exhibit any kind of collective behaviour as 
seen in these heavy-ion collisions. An indication for this might be the observation of long range, near- 
side angular correlations (ridge) in pp collisions at ^/s = 0.9, 2.36 and 7 TeV with charged particle 
multiplicities above four times the mean multiplicity 1221 f23l . Since quarkonium yields in heavy-ion 
reactions aie expected to be modified relative to minimum bias pp collisions Il24ll25l l26 l. one might ask 
whether their production rates in high-multiplicity pp colUsions are already exhibiting any effect like 
J/ V/^ suppression. 

In this Letter, we report the first measurement of relative J/y production yields {dNj/xf, / dy) / {dNj/xf,/dy) 
at mid-rapidity {\y\ < 0.9) and at forward rapidity (2.5 < j < 4) as a function of the relative charged 
particle multiplicity density (dA'ch/dT])/(dA'^ch/dT7) as determined in |t7| < 1 for pp collisions at ^/s = 
7 TeV at the LHC. 

The data discussed here are measured in two complementary parts of the experimental setup of AL- 
ICE ETl : the central barrel (|t7L< 0.9) for the J/y detection in the di-electron channel and the muon 
spectrometer (—4 < Tj < — 2.5)ljfor J/i/a — ix^fx^ measurements. 

The central barrel provides momentum measurement for charged particles with > 100 MeV/c and 

'in the ALICE reference frame the muon spectrometer is located at negative z positions and thus negative 
(pseudo-)rapidities. Since pp collisions are symmetric relative to 3; = 0, we have dropped the minus sign when rapidities 
are quoted. 
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particle identification up to pi « 10 GeV/c. Its detectors are all located inside a large solenoidal magnet 
with a field strength of 0.5 T. Used in this analysis are the Inner Tracking System (ITS) and the Time 
Projection Chamber (TPC). The ITS Il28l consists of six layers of silicon detectors surrounding the beam 
pipe at radial positions between 3.9 cm and 43.0 cm. Silicon Pixel Detectors (SPD) aie used for its 
innermost two layers and allow a precise reconstruction of the interaction vertex. The TPC |29 | is a 
large cylindrical drift volume covering the range along the beam axis relative to the Interaction Point (IP) 
between —250 < z < 250 cm and extending in the radial direction from 85 cm to 247 cm. It is the main 
tracking device in the central barrel and is also used for particle identification via a measurement of the 
specific ionization {AE /dx) in the detector gas with a resolution of about 5% ETl . 

The muon spectrometer consists of a frontal absorber followed by a 3 T-m dipole magnet, coupled to 
tracking and triggering detectors. Muons are filtered by the 10 interaction length (Ai) thick absorber 
placed between 0.9 m and 5.0 m from the nominal position of the IP along the beam axis. Muon tracking 
is performed by five tracking stations, positioned between 5.2 m and 14.4 m from the IP, each consisting 
of two planes of cathode pad chambers. The muon triggering system consists of two stations positioned 
at 16.1 m and 17.1 m from the IP, each equipped with two planes of resistive plate chambers. It is located 
downstream of a 1.2 m thick iron wall (7.2 X\) which absorbs hadrons penetrating the frontal absorber, 
secondary hadrons escaping the absorber material, and low-momentum muons {p <A GeV/c). A conical 
absorber surrounding the beam pipe provides protection against secondary particles throughout the full 
length of the muon spectrometer. 

Two VZERO detectors are used for triggering on inelastic pp interactions and for the rejection of beam- 
gas events. They consist of scintillator an^ays and are positioned at z = —90 cm and z = +340 cm, 
covering the pseudorapidity ranges —3.7 < T] < — 1.7 and 2.8 < Tj < 5.1. The minimum bias (MB) pp 
trigger uses the information of the VZERO detectors and the SPD. It is defined as the logical OR between 
two conditions: (i) a signal in at least one of the two VZERO detectors has been measured; (ii) at least 
one readout chip in the SPD fires. It has to be in coincidence with the arrival of proton bunches from 
both sides of the interaction region. The efficiency of the MB trigger to record inelastic collisions was 
evaluated by Monte Carlo studies and is 86.4% |[30l . For the di-muon analysis, a more restrictive trigger 
is used (/I -MB). It requires the detection of at least one muon above a threshold of pj"^ > 0.5 GeV/c in 
the muon trigger chambers in addition to the MB trigger requirement. 

The results presented in this Letter are obtained by analyzing pp collisions at = 7 TeV recorded in 
2010. For the J/i/a measurement in the di-electron (di-muon) channel a sample of 3.0 x 10^ minimum 
bias events (6.75 x 10^ pL-MR triggered events) is analysed, corresponding to an integrated luminosity 
of 4.5 nb^' (7.7 nb^'). The di-electron sample is divided into four sepai^ate data sets with slightly 
different running conditions. The relative normalization between the number of /i-MB and minimum 
bias triggers needed to extract the integrated luminosity in the di-muon case is calculated using the 
ratio of the number of corresponding single muons with pt > 1 GeV/c. The luminosity at the ALICE 
interaction point was kept between 0.6 and 1.2 -10^^ cm^^ s^' for all the data used in this analysis. This 
ensures a collision pile-up rate of 4%, or lower, in each bunch crossing. In the case of the di-muon 
analysis the interaction vertex is reconstructed using tracklets which are defined as combinations of two 
hits in the SPD layers of the ITS, one hit in the inner layer and one in the outer. Since for MB trigger 
used in the di-electron analysis the full information of the central barrel detectors is available (/i-MB 
triggered events only include SPD information), tracks measured with ITS and TPC are used in this case 
to locate the interaction vertex. This results in a resolution in z direction of ~ 600 /N^^^ /xm, where A^trk 
is the multiplicity measured via SPD tracklets. For the vertices reconstructed using SPD tracklets only, 
this resolution is worse by 35% for high (A/^trk = 40) and 50% for low (A^trk = 10) multiplicities. Events 
that do not have an interaction vertex within |zvtx| < 10 cm are rejected, where Zvtx is the reconstructed 
z position of the vertex. The rms of the vertex distributions along z is for all running conditions below 
6.6 cm. 
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Fig. 1: The distribution of the relative charged particle density {dNch/drj)/ (d/Vch/dfj) reconstructed around mid- 
rapidity (|tj| < 1.0) after correction for SPD inefficiencies. The vertical lines indicate the boundaries of the multi- 
plicity intervals used in this analysis. 



Pile-up events are identified by the presence of two interaction vertices reconstructed with the SPD. 
They are rejected if the distance along the beam axis between the two vertices is larger than 0.8 cm, 
and if both vertices have at least three associated tracklets. This removes 48% of the pile-up events. 
In the remaining cases two events can be merged into a single one, thus yielding a biased multiplicity 
estimation. A simulation assuming a Gaussian distribution for the vertex z position results in a probability 
for the occurrence of two vertices closer than 0.8 cm of 7%. Combined with the pile-up rate of 4%, this 
gives an overall probability that two piled-up events are merged into a single event of 0.3%, which is 
a negligible contribution in the multiplicity ranges considered here. 

The charged particle density dNch/'i'n is calculated using the number of tracklets Ntrk reconstructed from 
hits in the SPD detector, because the SPD is the only central baiTcl detector that is read out for all 
of the jLi-MB trigger. The tracklets are required to point to the reconstructed interaction vertex within 
±1 cm in radial and ±3 cm in z direction |[3Tl[32ll . Using simulated events, it is verified that A'trk is pro- 
portional to d/Vch/dT]. For a good geometrical coverage, only tracklets within |t7| < 1 from events with 
|Zvtx| < 10 cm are considered. Since the pseudorapidity coverage of the SPD changes with the interaction 
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Table 1: The boundaries of the used charged particle multiplicity intervals as defined via the number of SPD 
tracklets A^trk' the corresponding charged particles density ranges and mean values (dA^|,h/dJ?(bin)), as well as the 
number of analyzed minimum bias triggered events in the di-electron (A^evt'^ ) and the di-muon channel (A^^^ ). 
In the latter case this is the equivalent number of events, derived from the number of ji-MB triggered events. 

Multiplicity interval Mrk interval dA/ch/d^] range (dA/^ch/dT](bin)) N^^^";' x 10'' A'^fqlevt. x 



1 [1,8] 0.7-5.9 2.7 164.6 262.0 

2 [9,13] 5.9-9.2 7.1 51.1 79.5 

3 [14,19] 9.2- 13.2 10.7 35.7 55.4 

4 [20,30] 13.2-20.4 15.8 28.5 44.4 

5 [31,49] 20.4-32.9 24.1 9.7 15.3 



vertex z position and also with time, due to the varying number of dead channels, a correction is deter- 
mined event-by-event from measured data as a function of Zvtx and for each analysed data set separately. 
Figure [T] shows the resulting distribution of the relative charged particle density (dA'^ch/dTj)/(dA'^ch/dTj), 
where (dNch/dTl) = 6.01 it 0.01(stat.)^Q j2(syst.) as measured for inelastic pp collisions with at least 
one charged particle in |t]| < 1 [32]. The use of relative quantities was chosen in order to facilitate the 
comparison to other experiments and to theoretical models, as well as to minimize systematic uncertain- 
ties. The definition of the charged particle multiplicity intervals used in this analysis is given in Table [T] 
together with the corresponding mean values of dNch/drj. The present statistics allows one to cover 
charged particle densities up to four times the minimum bias value. 

For the I/i/A measurement in the di-electron channel tracks are selected by requiring a minimum pt of 
1 GeV/c, a pseudorapidity range of |t]| < 0.9, at least 70 out of possible 159 points reconstructed in 
the TPC and an upper limit on the l^-d-f- from the momentum fit of 2.0. Furthermore, tracks that 
are not pointing back to the primary interaction vertex within 1.0 cm in the transverse plane and within 
3.0 cm in z direction are discarded. To further reduce the background from conversion electrons a hit 
in at least one of the four innermost ITS layers is also required. Particle identification is performed 
by measuring the specific ionization AE/Ax in the TPC. All tracks within ib3a around the expected 
AE /Ax signal for electrons and at the same time outside ib3a (±3. 5a) around the expectation for protons 
(pions) are accepted as electron and positron candidates. e+ and e^ candidates that form a pair with 
any other candidate with an invariant mass below 0.1 GeV/c^ are discarded to reduce the amount of 
electrons coming from / conversions or Dalitz decays as well as their contribution to the combinatorial 
background in the di-electron invariant mass spectrum. 

The invariant mass distributions of the e+e^ pairs are recorded in intervals of the charged particle multi- 
plicity as measured using the SPD tracklets. As an example, the lowest and highest multiplicity intervals 
are shown in the two left panels of Fig. [2] The combinatorial background in each multiplicity interval is 
well described by the track rotation method, which consists in rotating one of the tracks of a e+e^ pair 
measured in a given event around the z axis by a random -angle in order to remove any conelations. 
After subtracting the background, the uncorrected I/i// yields are obtained by integrating the distribution 
in the mass range 2.92 - 3.16 GeV/c^. Alternatively, the combinatorial background is estimated by like- 
sign distributions, N^^ +N . These are scaled to match the integral of the opposite-sign distributions 
in the mass range above the J/ i/A signal (3.2 < m-mv < 4.9 GeV/c^) in order to also account for correlated 
background contributions, which mainly originates from semi-leptonic charm decays. Both methods 
provide a good description of the combinatorial background and their comparison is used to evaluate the 
systematic uncertainty on the 1/ 1// signal. 

For the I/i/a analysis in the di-muon channel muon candidates are selected by requiring that at least one 
of the two muon candidates matches a trigger track reconstructed from at least three hits in the trigger 
chambers. This efficiently rejects hadrons produced in the frontal absorber and then absorbed by the iron 
wall positioned in front of the trigger chambers. Furthermore, a cut /?abs > 17.5 cm is applied, where 
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Fig. 2: Opposite sign invariant mass spectra of the selected electron [(a) + (c)] and muon [(b) + (d)] pairs (filled 
symbols) for the lowest [(a) + (b)] and highest [(c) + (d)] multiplicity intervals. Also shown are the estimates of 
the combinatorial background which are based on a fit to the fl^fl^ pair distributions (solid line), and on like-sign 
pairs (open circles), as well as track rotation (open squares), in the e+e^ case. The number of events quoted in the 
figures refer to the corresponding minimum bias triggered events. 



Rabs is the radial coordinate of the track at the end of the frontal absorber (z = —5.03 m). Such a cut 
removes muons produced at small angles that have crossed a significant fraction of the thick beam shield. 
Finally, a cut on the pair rapidity (2.5 < j < 4) is appUed to reject events very close to the edge of the 
spectrometer acceptance. 

The number of J/i/a in each multiplicity interval is obtained by fitting the corresponding di-muon in- 
variant mass distribution in the range 2 < m-mv < 5 GeV/c^. The line shapes of the J/xjf and V^(2S) are 
parametrized using Crystal Ball functions f33l, while the underlying continuum is fitted with the sum 
of two exponential functions. The parameters of the Crystal Ball functions are adjusted to the mass dis- 
tribution of a Monte Carlo signal sample, obtained by generating J/y and i//^(2S) events with realistic 
phase space distributions [8 ]. Apart from the J/y/ and i//^(2S) signal normalization, only the position of 
the J/i/A mass pole, as well as its width, are kept as free parameters in the fit. Due to the small statistics, 
the V'^(2S) mass and width are tied to those of the J/i/a, imposing the mass difference between the two 
states to be equal to the one given by the Particle Data Group (PDG) [34 ], and the ratio of the resonance 
widths to be equal to the one obtained by analyzing reconstructed Monte Carlo events. Details on the 
fit technique can be found in [8|. The two right panels of Fig. |2] show the measured di-muon invariant 
mass distributions together with the results of the fit procedure for the lowest and highest multiplicity 
intervals. 

The results are presented as the ratios of the J/y yield in a given multiplicity interval relative to the 
minimum bias yield. By performing simulation studies in intervals of dNch/drj it was verified that the 
geometrical acceptances, as well as the reconstruction efficiencies, do not depend on dNch/'^'H in 
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Fig. 3: J/y/ yield cWj/^/dy as a function of the charged particle multiplicity densities at mid-rapidity d^Vj-h/di?- 
Both values are normalized by the corresponding value for minimum bias pp collisions ((dA'j/^/dy), {dNch/drj))- 
Shown are measurements at forward rapidities ji^pL^ , 2.5 < y < 4) and at mid-rapidity (J/\// e+e^, 

\y\ < 0.9). The error bars represent the statistical uncertainty on the J/xj/ yields, while the quadratic sum of the 
point-by-point systematic uncertainties on the J/xj/ yield as well as on d^Vj-i^/dTj is depicted as boxes. 



range under consideration here (dA^ch/dT] < 32.9). Therefore, these corrections and their correspond- 
ing systematic uncertainties cancel in the ratio {dNj/xf,/dy)/{dN]/xf,/dy) and only the uncorrected signal 
counts have to be divided. The number of events used for the normalization of {dN]/^,/dy) is corrected 
for the fraction of inelastic events not seen by the MB trigger condition. After applying acceptance and 
efficiency corrections these values are in agreement with those that can be obtained from the numbers 
quoted in [8]: (dNy^/dy) = (8.2 ± 0.8(stat.) ± 1.2(syst.)) x 10"^ for J/y ^ e+e" in \y\ < 0.9, and 
{dNy^/dy) = (5.8±0.2(stat.)±0.6(syst.)) x 10"^ for J/i/a ^ in 2.5 <y<4. In the case of the 
3/Y yields measured in a given multiplicity interval, no trigger-related correction is needed, since the 
trigger efficiency is 100% for A^trk > 1- 

The systematic uncertainties are estimated as follows. In case of the di-electron analysis, the absolute 
differences between the resulting {dN]/xf,/dy) / {dNj/xf,/dy) values obtained by using the like-sign and the 
track rotation methods define the uncertainty due to the background subtraction. It is found to vary 
between 2% and 12% for the different multiplicity intervals. For the di-muon analysis this uncertainty 
is evaluated by varying the functional form of the background description (polynomial instead of sum of 
two exponential). It depends on the signal to background ratio and varies between 3% and 4%. Since 
for the muon measurement it is not possible to associate a measured track to the interaction vertex, an 
additional systematic uncertainty arises from pile-up events. Among the vertices inside these events 
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Fig. 4: Relative J/y/ yield dN^^/dy as a function of relative charged particle multiplicity densities around mid- 
rapidity dA^ch/di? as calculated with PYTHIA 6.4 in the Perugia 201 1 tune ll35l[36ll . Shown are results for directly 
produced J/ i/A in hard scatterings via the NRQCD framework at forward rapidities (2.5 < y < 4) and at mid-rapidity 
(bl<0.9). 



always the one with the largest number of associated tracks is chosen as main vertex. Therefore, events 
with very low multiplicities are more likely to have a wrong assignment and thus this uncertainty is largest 
in the first multiplicity interval (6%), while it is 3% in the others. Possible changes of the spectra 
with event multiplicity can introduce a dA/^ch/drj dependence of the acceptance and efficiency correction, 
thus resulting in an additional systematic uncertainty. This is estimated by varying the {p^) of the J/i/a 
spectrum that is used as input to the determination of the corrections via simulation between 2.6 and 
3.2 GeV/c. A systematic effect of 1.5% (3.5%) is found for the di-electron (di-muon) analysis. The total 
systematic error on {ANn^, / Ay) / {ANn^, / Ay) is given by the quadratic sum of the separated contributions 
and amounts to 2.5 -12% depending on the multiplicity interval for the di-electron result. In the case 
of the di-muon analysis it varies between 8% in the first and 6% in the last multiplicity interval. An 
additional global uncertainty of 1.5% on the normalization of {ANii^,/Ay) is introduced by the correction 
of the trigger inefficiency for all inelastic collisions. 

The systematic uncertainties on {ALNch/Ar\)/{&Nc]x/Ar]) are due to deviations from a linear dependence 
of ANch/Ar\ on Nti\^ and variations in the Ni^c^ distributions which remain after the correction procedure. 
The latter are caused by changes in the SPD acceptance for the different data taking periods. The first 
contribution is estimated to be 5%, while the second is ~ 2%, as determined by Monte Carlo studies. In 
addition, the systematic uncertainty of the [ANch/Ax]) measurement (^2 0%) 1321 is also included. 

Figure [3] shows the relative J/i/a yields measured at forward and at mid-rapidity as a function of the 
relative charged particle density around mid-rapidity. An approximately linear increase of the relative 
J/i/A yield {dNy^i/dy)/ {dNn-^/dy) with (dA/^ch/dT7)/(dA'^ch/dT7) is observed in both rapidity ranges. The 
enhancement relative to minimum bias J/i//^ yield is a factor of approximately 5 at 2.5 < j < 4 (8 at 
\y\ < 0.9) for events with four times the minimum bias charged particle multiplicity density. 

A possible interpretation of the observed correlation of the J/ yield and the charged particle multiplicity 
might be that J/i// production is always accompanied by a strong hadronic activity, thus biasing the 
dA'^ch/dTj distributions to higher values. Since this correlation extends over the three units of rapidity 
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between the mid-rapidity dA'^^h /d^] and the forward rapidity J/ i/a measurement, it would have far reaching 
consequences on any model trying to describe J/ i/a production in pp collisions. 

In order to illustrate that the observed behaviour cannot be understood by a simple 2 — )• 2 hard partonic 
scattering scenario, a prediction by PYTHIA 6.4.25 in the Perugia 201 1 tune [35 , 36] is shown in Fig.|4]as 
an example. Only J/ directly produced in hard scatterings via the NRQCD framework |[37l (MSEL=63) 
are considered, whereas J/ 1//^ resulting from the cluster formation processes are ignored. A J/i//^ cluster 
is a string formed by a cc pair produced via parton shower evolution which has an invariant mass that 
is too low for the standard Lund string fragmentation procedure and thus does not correspond to a well 
defined hard scattering process. The calculation shown in Fig. |4] is thus the ratio of the multiplicity 
distributions generated for minimum bias events and events containing J/i/a from hard scatterings. It 
exhibits a decrease of the J/i/a multiplicity with respect to the event multiplicity, which indicates that hard 
J/ production, as modeled by PYTHIA 6.4.25, is not accompanied by an increase of the total hadronic 
activity. Further studies with other models such as PYTHIA 8 138)1 and Cascade |[391 are needed. It 
should be pointed out that the measurement also includes J/i/a from the decay of beauty hadrons, which 
is not part of the shown PYTHIA result. The fraction of J/y from feed down can change with the 
event multiplicity and can therefore contribute to the observed multiplicity dependence. However, since 
this contribution is on the order of 10% E Ul IHl it might be only a small contribution to the observed 
differences between model and data. 

On the other hand, the increase of the J/ i/A production with event multiplicity, as reported here, might be 
due to MPI. In this scenario the multiplicity of charged particles is a direct measurement of the number 
of partonic interactions in the pp events. If the effect of MPI extends into the regime of hard processes, 
also the J/ 1//^ yield should scale with the number of partonic collisions and the observed correlation will 
result. It has even been conjectured in HOTl that the increase of the J/yA yield with dA^ch/dTj and the 
ridge phenomenon observed in high-multiplicity pp collisions f23] could be related. They might both be 
caused by gluon density fluctuations and a special transverse distribution of the gluon densities inside 
the nucleon. The presence of these fluctuations could significantly increase the probability for MPI and 
thus cause the observed rise of the J/i/a yield. 

The multiplicity dependence measured here will allow a direct comparison of the J/y production in pp to 
the one observed in heavy-ion collisions. With a mean value of dA/^ch/dr] of 24.1, the highest multiplicity 
interval shown in Fig. [3j for instance, corresponds roughly to 45 - 50% centrality for Cu-Cu collisions 
at ^/s^ = 200 GeV |[2TI . In order to establish whether any evidence for a J/i/A suppression is observed 
already in pp, a proper normalization is needed. This could be provided by a measurement of open charm 
production in the same multiplicity bins. Corresponding studies are currently ongoing. 

In summary, yields are measured for the first time in pp collisions as a function of the charged 
particle multiplicity density dNch/drj. J/y mesons are detected at mid-rapidity {\y\ < 0.9) and forward 
rapidity (2.5 <y <4), while dNch/drj is determined at mid-rapidity (|t7| < 1). An approximately linear 
increase of the J/y yields with the charged particle multiplicity is observed. The increase is similar at 
forward and mid-rapidity, exhibiting an enhancement relative to minimum bias J/y yield by a factor 
of about 5 at 2.5 < y < 4 (8 at |y| < 0.9) for events with four times the minimum bias charged particle 
multiplicity. Our result might either indicate that J/y production in pp collisions is always connected 
with a strong hadronic activity, or that multi-parton interactions could also affect the harder momentum 
scales relevant for quarkonia production. Further studies of charged particle multiplicity dependence of 
J/VA, T, and open charm production, also as a function of pi, will shed more light on the nature of the 
observed effect. 
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